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Abstract 

The orbital period of the binary pulsar PSR 1913+16 has been observed 



| to decrease at the rate of 2.40 x 10 12 s/s which agrees with the prediction 

Q\ ' of the quadropole formula for gravitational radiation to within one percent. 

| The decrease in orbital period may also occur by radiation of other massless 

^ ' particles like scalars and pseudoscalar Nambu-Goldstone bosons. Assuming 

o ' 

CN ' that this energy loss is less than one percent of the gravitational radiation, 

o 



we can establish bounds on couplings of these particles to nucleons. For a 
scalar nucleon coupling of the form g s <f) ip if) we find that g s < 3 x 10~ 19 . 
From the radiation loss of massless Goldstone bosons we establish the upper 



bound 0/f< 7.5 x 10 GeV on the QCD vacuum angle 9 and the scale 
/ at which the baryon number symmetry can be broken spontaneously. 
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The Hulse- Taylor (H-T) binary consisting of the pulsar PSR 1913+16 orbiting around 
an unseen companion star provides firm evidence for the existence of gravitational waves 
]1||. The observed loss of orbital period agrees with the prediction from the quadropole 
formula of gravitation radiation to within one percent. In this letter we compute the 
orbital energy loss due to radiation of other massless particles like scalars and pseudoscalar 
Nambu-Goldstone bosons. Massless scalars which couple to nucleons arise in scalar-tensor 
theories of gravity |§, as dilatons in theories with spontaneously broken conformal symmetry 
|H and in string theories || . Assuming a generic scalar- nucleon coupling C s = g s <p s ip if) we 
find that the radiation of <fi s particles from the H-T binary is less than 1% of the gravitational 
radiation loss if g s < 3x 10~ 19 . This gives an upper bound a# = g 2 s / 'AnGrn^ < 1 on the ratio 
of a long range scalar mediated fifth force to the gravitational force between two nucleons. 
This bound is not as stringent as the bounds obtained from terrestrial fifth force search 
experiments which give as in the range of 10~ 3 — 10 -6 . 

In theories with a spontaneously broken global symmetry like the baryon number or 
the lepton number we have massless Nambu-Goldstone bosons (NGB) which have a generic 
coupling C p = (m/ f)(j) p ipi^^ip where m is the fermion mass and / is the scale of the global 
symmetry breaking ||. The pseudoscalar field of a macroscopic source adds up coherently 
only if the spins of the constituents are polarised. It was observed by Chang, Mohapatra and 
Nussinov ]l(J and Barbieri et al that the CP violating operator 9GG in the QCD sector 
induces a coupling C = (9 / f){m u m d / (m u + m d ))(j) p ^ip between the NGB cf> p and nucleons. 
This coupling give rise to a 1/r type long range force and the NGB field of the constituent 
nucleons a macroscopic test body add up coherently even when their spins are randomly 
aligned. From the constraints on the energy carried away by the radiation of NGBs from the 
H-T binary we obtain the upper bound (8 / f) < 7.5 x 10~ 16 GeV -1 . This can be compared 
with the separately established bounds 9 < 1CT 9 (from the measurement of the neutron 



electric dipole moment fi2| , |13|| ) and / > 10 GeV (from the cooling rate of helium stars 
The rate of energy loss by scalar particle emission is oc f2 4 (where Q is the orbital 
frequency). Observations of binary systems with faster orbital frequencies can be used to 



put more stringent bounds on couplings considered in this paper. 

Finally we compute the energy loss by the radiation of neutrino pairs from the con- 
stituent neutrons of the H-T binary. We find that for the neutral current coupling 
x)w, the energy radiated by neutrino pair emission is suppressed by 
the phase factor and is negligibly smaller than the gravitational radiation. Therefore if 
experimentally one observes a discrepancy between the observed period loss of the binary 
orbit and the prediction from gravitational radiation formula, it would be a signal of new 
kind of massless particle radiation and signal of new physics beyond the standard model. 



Gravitational radiation : We use the Feynman rules of linearised quantum gravity |lq -|T8 



to compute the gravitational radiation from the Hulse- Taylor binary system. Assuming a 



universal graviton-matter coupling Cj = v8ttG h^ u T^ v we find that the energy loss by 
gravitational brehmstrahlung at tree level agrees with the Peter-Mathew expression for 
classical gravitational wave radiation from binary system. 

The effective Lagrangian for the graviton matter interaction is 

C = -h^nh llu + K~h^ (1) 

where the graviton field h^ u is a perturbation of the metric g^ v = rj^ + nh^ to the first 
order in k = y/8irG. We have made the harmonic gauge choice, dfji^ \ = (l/2)d u h^ and 
defined h^ v = h^ u — \rj^ u h x \ The universal coupling of W with matter at the tree level is 
a reflection of the equivalence principle of classical gravity. The rate of graviton emission is 
given by 

2 d 3 k 1 

dT = ^ £ \T» v {k>) eft (k)\ 2 2?r S(u - u/) ^ - (2) 

where T^ik') and ^\{k) are the Fourier transforms of T^ v (x) and W(x) and A denotes 
the polarisation of the emitted gravitons which is summed over. Using the harmonic gauge 
condition /c M (e^ — ~ r]^e x \) = 0, we can set e M o = 0, e l , = 0, k l eij = we have e^ v = 
equation (fj) can be written as 

dT = \ (27p ? e W* (fc) e W (A ° T: ^' ] Tlm{u),) ^ 6{U " U>) dU dQk (3) 



Using the relation 



57T 



A=l 



.-, ydudjm + 5i m 8jl) — — %5j m 



(4) 



We have the graviton emission rate 



dr 



57T 



T^TUu') - -\T 1 t (u')\ 2 ]uj S(uj - u') du 

o J 



(5) 



and the rate of energy loss by the graviton radiation is 



dE 
~dt 



5tt 



T^T^J) - -\T * ,(u/)| 2 5(u - u') du 



(6) 



For a binary system of stars with masses mi, m 2 in an eliptical orbit around the center of 
mass, the stress tensor is given by 



T» v {x) = M8\x-y{t))UJJ v 



(7) 



where M = m 1 m 2 /(m 1 + m 2 ) is the reduced mass and = 7(1, x, y, z) is the four 
velocity of the reduced mass in the elliptical orbit. Assuming an orbit in the x — y plane 
with coordinates y(t) = (d cos8, d sinO), the relative distance d and the angular velocity 9 
are given by 

a(l -e 2 ) 



d 



1 + e cosO 



G(m\ + m 2 ) a (1 — e 2 ) 
I 2 



1/2 



(8) 



(9) 



where a is the semimajor axis and e the eccentricity of the eliptical orbit. Substitution @ 
and (H) in (0) and taking the fourier transform we obtain the stress tensor Tij(u' = nfl) 
for the Kepler orbit [2] in terms of the n harmonics of the fundamental frequency Q = 
(G(mi + m 2 ) a" 3 ) 1/2 , 



,/2 



T xx (uj') = M a 2 ^[J n _ 2 (ne) - 2e J n _i(ne) + 2e J n+1 (ne) - J n+2 {ne)] (10a) 



Tyy(uj') 



u- . , 2 2 



— M a 2 -[J n _ 2 (ne) 
4 n 



- 2e J n _i(ne) + - J n (ne) + 2e J n+ i(ne) - J n+2 (ne)] (10b) 



J* 2 

T xy {u') = — M a 2 — (1 — e 2 ) 1/2 [J n _ 2 (ne) - 2 J n {ne) + J n+2 {ne)} (10c) 
From the above equation we obtain 

T,>')W) - 1 ^')| 2 ] = a 4 a/ 4 <?(n, e ) (11) 

where 

e) = 1 {[J n _ 2 (ne) - 2e J n _i(ne) + 2e J n+ i(ne) - J n+2 (ne)] 2 

+ (1 - e 2 )[J r ^ 2 {ne) - 2 J n (ne) + J n+2 (ne)} 2 + -^[J n {ne)f } (12) 

Substituting ([□]) into @ we have the rate of energy loss of a binary system in an elliptical 
orbit as sum of radiation in the n harmonics of the fundamental frequency Q, 

^ = ^ £(« «)" ■ M 2 a 4 (n fi) 4 ^(n, e) (13) 



Using the relation ™ 6 e) = (1 - e 2 )~ 7/2 (l + ||e 2 + §|e 4 J [from [g]] and substituting 
the values of Q = 0.2251 x 10~ 3 sec" 1 , m x = 1.42M , m 2 = 1.4M , a = 3.0813815 lsec, 
e = 0.617127 for the parameters of the H-T binary, we find that the energy loss is 



— = — ■ G ■ fi 6 ■ ( mim2 ) 2 a 4 (1 - e 2 )^ 2 ( 1 + ^e 2 + ^ 
eft 5 Vmi+m 2 / v ; V 24 96 



= 3.2 x 10 33 erg/sec (14) 

The time period of the elliptical orbit depends upon the energy E, so energy loss leads to a 
change of the time period of the orbit at the rate 

^ = -67rGHVVima)~Vi + m 2 )" 1/2 a 5 / 2 ^ 

192 ^5/3 Q5 /3 ™i™2 _ e2) „ 7/2 / + 73 2 + 37 A g) 



5 (mi+m 2 ) 1 /3 v y V 24 96 

which is the same expression as obtained by Peter and Mathews |2| for the classical grav- 
itational radiation from binary systems. We find therefore that the tree level quantum 
gravity calculation agrees with the clasical results. For the Hulse - Taylor binary, the 
expression (|I5|) yields the orbital period decceleration due to the gravitational radiation 



Pb = —2.403 ± 0.002 x 10~ 12 which agrees with the observed value from the Hulse - Taylor 
binary |lj Pt,(observed) = —2.40 ± 0.09 x 10~ 12 to with in 1%. Energy loss by emission 
of other massless particles should be within 1% of energy loss by gravitational radiation 
and that can be used to put bounds on the couplings of the various massless scalar and 
pseudoscalar particles which arise in particle physics models. 

Massless scalar radiation: We assume a coupling between massless scalar fields <p s and the 
baryons of the form 

£ s = g s 4>s i> ip (16) 

which for a macroscopic baryon source can be written as 

C s = g s 4> s n(x) (17) 

where n(x) is the baryon number density. A neutron star with radius ~ 10 km can be 
regarded as a point source since the Compton wavelength of the radiation ~ Q^ 1 = 10 9 km 
is much larger than the dimension of the source. The baryon number density n(x) for the 
binary stars (denoted by a = 1, 2) may be written as 

n(x) = N » S 3 (x-x a (t)) (18) 

a=l,2 

where N a ~ 10 57 is the total number of baryons in the neutron star and x a (t) represents the 
Keplerian orbit of the binary stars. For the coupling flT7| ) and the source (18) the rate of 
scalar particles emitted from the neutron star in orbit with frequency Q is 

dT = |nM| 2 (2vr) 5{u> - u') (19) 

the rate of energy loss by massless scalar radiation is 

(20) 

where n(u) is the fourier expansion of the source density (18) 

;)= — / ^ e - iwt V N a 5 3 (x-x a (t))d 3 x dt (21) 

27r J a=l,2 



n[uj) 



with u = nQ. Going over to the cm. coordinates r 
X\ = m 2 rj (mi + m 2 ), x 2 = — m\rj (mi + m 2 ) we have 



(x, y) by substituting 



n(u) = (JVi + N 2 ) 5(cu) + (— - —) M ( ik x x{u) + ik y y(u)) + 0{k, ff (22) 

\ 777,1 777,2/ 

where (x(u),y(uj)) are the fourier components of the Kepler orbit of the reduced mass in the 
cm. frame given by (8)-(9). The first term in n(u) is a delta function which has vanishing 
contribution to (20). The leading non-zero contribution comes from the second term in (22). 
Substituting || 



2 a 2i a y/1 - e 
x{u) = — J n {ne), y{u) = J n {ne) 

71/ Tl/G 



(23) 



in (22) we obtain the expression for |n(c<j)| 2 given by 



\n{u)f 



Nt N 2 



M 



1 



a 2 Q 2 [J'l(ne) + ^ ' * 



J 2 n (ne)] 



(24) 



,777i 777,2. 

where we have used the dispersion relation /c 2 = k 2 y = |(nf2) 2 . Substituting (24) in (20) 
we have the rate of energy loss by massless scalars 



dE 



2 

3^ 



M g s 

mi m 2 , 



^ E ^ [J nine) + 



2, , , (1-e 2 ; 



J>e)\ 



(25) 



The mode sum can be carried out using the Bessels function series formulas given in reference 
§ to give En n 2 [J'l{ne) + (1 - e 2 )e~ 2 J 2 (ne)} = (1/4) (2 + e 2 )(l - e 2 )" 5 / 2 . The energy 
loss (25) in terms of the orbital parameters Q, a and e is given by 



dE 
~dT 



1 

3^ 



m n 2 

M g s 

mi 777,2- 



1 - e 2 ) 5 /2 



(26) 



Since N a m n = m a — e a where e 
neutron mass, the factor ( — — — 

' \ m l m 2 



Gm 2 n 



„ - is the gravitational binding energy and m n the 

G - if ) • For the H_T binar y o.o2M Q 

and R a ~ 10 km, therefore f ^ - ~ 3 x 10~ 3 GeV" 1 . 

u ' \ mi ni2 / 



For the H-T binary the rate of energy loss terms out to be 



dE 



= g 2 s x 9.62 x 10 67 ergs/sec 

Lblj 



(27) 



Assuming that this is less than 1% of the gravitational energy loss i.e. dE/dt < 10 31 
ergs/sec, we obtain gives an upper bound on scalar nucleon coupling 

g s < 3 x KT 19 . (28) 

Exchange of massless scalar between two nucleons gives rise to spin independent fifth force 
with the static potential V ss (f) = — g^/Airr, which shows that ratio of the fifth force to the 

2 

gravitational force between two nucleons is as = i7r Q m 2 < 1- This is less stringent than 
the bound as ~ 1CT 3 — 1CT 6 obtained from terrestrial fifth force search experiments [|. 
Nambu-Goldstone boson radiation: Massless pseudoscalar particles arise as Nambu- 
Goldstone bosons (NGB) when some global symmetry is broken spontaneously ||T9|-f2lH . We 
consider the coupling of massless NGBs <p p to baryons, 

I — 77? — 

c p = - (d,<p p )^Yi 5 ^ = j M&J^) (29) 

(where the last equality holds for fermions on mass shell) which arise in many particle 
physics models where a global symmetry is broken spontaneously at some scale / giving rise 
to Nambu-Goldstone bosons. This coupling gives rise to a spin dependent long-range force 
p2|] V(r) ~ (l//) 2 (l/r 3 )(o"i • 02 — 3(<7i • r)(°~2 • f)). The pseudoscalar coupling in the first 
order in g p is spin dependent and the field of a macroscopic body with randomly oriented 
constituents averages to zero. The radiation of NGB by the N constituent particles of the 
macroscopic system will be in coherent which means that the energy radiated by iV particles 
will be N times the single particle energy loss. This is different from scalar radiation which 
being coherent scales as N 2 times the single particle radiation. 

It was observed by Chang, Mohapatra and Nussinov Jl0| and Barbieri et al [|11|] that if 



there is CP violation in the theory then <p p can have both pseudoscalar coupling as in (|29| 



as well as scalar coupling of the form (|T6|). For example the CP violating QCD term 6GG 



will induce an interaction between (j) p and nucleons of the form 



Of m u md , 

£ P = 7 (30) 



/ \m u + m d 



This scalar coupling gives rise to a long range (V(r) ~ 1/r) potential which is spin inde- 
pendent so that cf) p field outside a macroscopic object adds up coherently. The form of the 
interaction term ([30]) is the same as in equation ( |16D so the same bound (28) holds for the 
dimensionless coupling 

/ m u m d \ < 3 x 1Q _ 19 
/ Vm u + m d y 

which means that for m u = 5 MeV, = 9 MeV we have 9/ f < 7.5 x 10~ 16 GeV -1 . This can 
be compared with the separate bounds 9 < 10~ 9 (obtained from neutron edm [|l^,|n|) and 



/ > 10 8 GeV (obtained from the cooling rate of helium stars The bound fl3T|) holds 



for models where the mass of the Nambu-Goldstone boson is smaller than the frequency 
Q ~ 10~ 19 eV of the binary orbit. In axion models p3l the axion acquires a mass due to 



QCD instantons which ranges between 10" 3 eV - 10" 6 eV, so our bound does not hold for 
these model. In majorons models |JT9| , f2"0"f where the lepton number is broken spontaneously, 



the majorons which arise as NGBs can remain massless and our bound would then hold. 
Neutrino radiation : The coupling of neutrinos to neutrons in the standard model is via the 
weak neutral current and is given by 

L v = V2G F ^-u Ll °u L (32) 

where n(x) is the number density of the neutrons which are the source of the neutrino field. 
The radiation of neutrinos from the H-T binary with n(x) = N 5 3 (x — x(t)) carries aways 
energy at the rate 

d£ „ (^i^j n 6 = 10- 43 ergs/sec. (33) 

The weak coupling of neutrinos to neutrons is much larger than the coupling of the gravitons 
to matter, however since the emission of neutrinos occurs in pairs, the phase space suppres- 
sion makes the energy radiated by neutrino emission negligible compared to graviational 
radiation. 

The period loss in the H-T system has been determined by measuring the time of pe- 
riastron over a period of almost 19 years. The accuracy of the measured value of period 

9 



loss increases quadratically with time. If in the course of observation one finds a significant 
discrepancy between the observed value of period loss and the prediction of the gravitational 
quadropole formula, it would be a compelling signal of physics beyond standard model. 
Note added: Astrophysical bounds from the Hulse- Taylor binary have been considered in 
references |M| and P5fl . In the energy loss in Brans-Dicke gravity is derived and in [25 
effect of scalar couplings on orbital parameters of the H-T binary is studied. 
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